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1. Introduction 

Laminar Hagen-Poiseuille and electro-osmotic flow is important to a variety of lab- 
on-a-chip applications and microfluidics E] and the rapid development of micro 
and nano fabrication techniques during the past decade has put even more emphasis 
on flow in channels with a variety of shapes depending on the fabrication technique 
in use. The list of examples includes rectangular channels obtained by hot embossing 
in polymer wafers, semi-circular channels in isotropically etched surfaces, triangular 
channels in KOH-etched silicon crystals, Gaussian-shaped channels in laser-ablated 
polymer films, and elliptic channels in stretched PDMS devices While general 
results for the shape-dependence of the hydraulic resistance in the case of a non- 
conducting fluid were reported recently there has, according to our knowledge, 
been no analogous detailed study of the shape-dependence of flow of electrolytes in 
the presence of a zeta potential which is a scenario of key importance to lab-on-a- 
chip applications involving biological liquids/samples in both microfluidic [HI El IE) and 
nanofluidic channels 131101 CIlinilTSlini. 

In this work we address the general steady-state flow problem in Fig. where 
pressure gradients and electro-osmosis (EO) are playing in concert |T5|. We consider 
a long, straight channel of length L having a constant cross section f2 of area A and 
boundary d£l of length V. The channel contains an incompressible electrolyte, which 
we for simplicity assume to be binary and symmetric, i.e., containing ions of charge 
+Ze and — Ze and equal diffusivities D. The electrolyte has the Debye screening 
length Ad, bulk conductivity cr , viscosity 77, permittivity e and at the boundary dfl 
it has a zeta potential The laminar, steady-state flow is driven by a linear pressure 
drop Ap and a linear voltage drop AV. With these definitions flow will be in the 
positive x direction. In the linear-response regime the corresponding volume flow rate 
Q and charge current I are related to the driving fields by 

AV h { G21 



Q 
I 



G 



G12 
G22 



(1) 



where, according to Onsager relations {Fg\, G is a symmetric, Gyi — G21, two-by-two 
conductance matrix. The upper diagonal element is the hydraulic conductance or 



p(0) = p Q + Ap 
V(0) = V + AV 




p{ l ) = Po 
V{L) = V 



Figure 1. A translation invariant channel of arbitrary cross section Q of area 
A containing an electrolyte driven by a pressure gradient —Ap/L and by electro- 
osmosis through the potential gradient —AV/L. The channel wall 90 has the 
electrical potential which induces a thin, charged Debye layer (dark gray) that 
surrounds the charge neutral bulk (light gray). 
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Gn 



(2) 



inverse hydraulic resistance given by 

A A 

arj L ' 

where a is the dimensionless geometrical correction factor, shown in Ref. |S] to be a 
linear function of the dimensionless compactness parameter C = V 2 / A. While there is 
no intrinsic length scale influencing Gn, the other elements of G depend on the Debye 
screening length A d ■ This length can be comparable to and even exceed the transverse 
dimensions in nano-channels [91 1101 111] , in which case the off-diagonal elements may 
depend strongly on the actual cross-sectional geometry. However, for thin Debye 
layers with a vanishing overlap all four matrix elements in G are independent of the 
details of the geometry. For a free electro-osmotic flow, a constant velocity field 
v co = (e(/r))AV/L is established throughout the channel, except for in the thin Debye 
layer of vanishing width. Hence Q = v co A and 



G12 — G21 — — - — , Xd *C — . 

rj L V 

From Ohm's law I = (a A/L)AV it follows that 



Lx92 — er — , 



A 

« t5- 



(3) 



(4) 



For strongly overlapping Debye layers we find, see Sections 13 . 2 . 21 and 14.21 that 
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2 sinh' 
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( Ze C 

\k B T 
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We emphasize that the above results are generally valid for symmetric electrolytes, 
even beyond the Debye-Hiickel approximation. In the Debye-Huckel limit ZeQ <^.ksT 
they also hold for asymmetric electrolytes. We also note that in the Debye-Hiickel 
limit the expressions agree fully with the corresponding limits for a circular cross 
section and the infinite parallel plate system, were explicit solutions exist in the Debye- 
Hiickel limit in terms of Bessel functions |171 IT%| and cosine hyperbolic functions ^21 j 
respectively. From the corresponding resistance matrix R — G _1 we get the hydraulic 
resistance 

r,L 



Rn 



a 



(7) 



where (3 = 
resistance. 



1 - p A 2 ' 

G12G21/ '(G11G22) is the Debye-layer correction factor to the hydraulic 
In the two limits we have 



2^2 



e 2 t 
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ZeC, 
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Xd« 15 
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(8) 



For C, going to zero /3 vanishes and we recover the result in Ref. [Hj for the hydraulic 
resistance. 
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2. Governing equations 

For the system illustrated in Fig. ^ an external pressure gradient Vp = — (Ap/L)e x 
and an external electrical field E = Ee x = (AV/L)e x is applied. There is full 
translation invariance along the x axis, from which it follows that the velocity field is 
of the form v(r) — v(r±)e x where r± = ye y + ze z . For the equilibrium potential and 
the corresponding charge density we have <fi cq (r) = 4> cq (r±) and p% q {r) = pt q (r±), 
respectively. We follow our related recent work ^] and use the Dirac bra-ket 
notation |2l)l 12 1| . where functions f{r±) in Q are written as |/) with inner products 
defined by the cross-section integral 

(f\g)= [ d? ± f<f ± )g{? ± ). (9) 

From the Navier-Stokes equation it follows that the velocity is governed by |221 |2~5| 

0= ^|l) +r/ vi|,) + ^|^ q ), (10) 

where = dy + &l is the 2D Laplacian. The equilibrium potential \4>eq) and the 
charge density \p% q ) are related by the Poisson equation 

Vl|&q)=~|pSq>- (11) 

The velocity \ v) is subject to a no-slip boundary condition on dVL while the equilibrium 
potential |</> e q) equals the zeta potential C on 9S1. Obviously, we also need a statistical 
model for the electrolyte, and in the subsequent sections we will use the Boltzmann 
model where the equilibrium potential |</> e q) is governed by the Poisson-Boltzmann 
equation. However, before turning to a specific model we will first derive general 
results which are independent of the description of the electrolyte. 

We first note that because Eq. (|10|) is linear we can decompose the velocity as 
ky = v p ) + 1 u C o) j where \v p J is the Hagen-Poiseuille pressure driven velocity governed 

by 

= ^|l)+r/Vi| Up ), (12) 
and l^eo) is the electro-osmotic velocity given by 

Ko} = -^(C|l}-|<£e q ». (13) 

The latter result is obtained by substituting Eq. (|llfl for |Pc q ) in Eq. IjlUI) . The upper 
diagonal element in G is given by G\\ — (l|u p )/Ap which may be parameterized 
according to Eq. (0). The upper off-diagonal element is given by G\2 = (l\v eo ) / AV 
and combined with the Onsager relation we get 

Gl2 = G 21 = - cq ) = -y-(C ~ &q), (14) 

where we have used that (l|l) = A and introduced the average potential <fi cq — 
<0 eq |l>/<l]l>. 

There are two contributions to the lower diagonal element G22] one from 
migration, G^ Ig = (l|<r)/£, and one from electro-osmotic convection of charge, 
G C 2 T = (pt q \v C o)/AV, so that 

G22 = G™ 8 + G C 2 T V = ^(lk) - ^<Pc q |C - eq ), (15) 
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where the electrical conductivity o~(r±) depends on the particular model for the 
electrolyte. For thin non-overlapping Debye layers we note that cq ~ so that 
Eq. I|14|) reduces to Eq. © and, similarly since the induced charge density is low, 
Eq. I|15fl reduces to Eq. Q . For strongly overlapping Debye layers the weak screening 
means that <p eq approaches £ so that the off-diagonal elements Gn = G21 and the 
G22 nv part of G22 vanish entirely. In the following we consider a particular model for 
the electrolyte and calculate the asymptotic suppression as a function of the Debye 
screening length A_d. 



3. Debye Hiickel approximation 

In the Debye-Huckel approximation the equilibrium potential \<t>eq) is governed by the 
linearized Poisson-Boltzmann equation |3] 

V lkcq) = T2-|0eq>, (16) 

where is the Debye screening length. The validity of this model will be discussed 
in more detail in Sec. 



3.1. Hilbert space formulation 

In order to solve Eqs. 1)11) [I . and (|16f) we will take advantage of the Hilbert space 
formulation |24|. often employed in quantum mechanics |21| . and recently employed 
by us on the problem of an accelerating Poiseuille flow \nj\ . The Hilbert space of real 
functions on £1 is defined by the inner product in Eq. @ and a complete, countable 
set (IV'n)} °f orthonormal basis functions, i.e., 

(ip m \ipn) = <W> (17) 

where S nm is the Kronecker delta. We choose the eigenf unctions {|VVi)} of the 
Hclmholtz equation (with a zero Dirichlet boundary condition on dfl) as our basis 
functions, 

- Vi|^ n ) = /4|^ n ), n= 1,2,3,.... (18) 

With this complete basis any function in the Hilbert space can be written as a linear 
combination of basis functions. In the following we write the fields as 



oo 



A 



) =^ a „|<0n}, (19) 

n=l 

oo 

'eq) =C|1)-5>«I^)> ( 20 ) 

oo 

eq)=E C «l^>- (21) 



Inserting Eqs. (|18fl and l|2(J|) into Eq. (IB) , and multiplying by (VVn|, yields 

K = C z . 2 , n = 1,2,3,... (22) 

1 + (KnXDr 
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Likewise, inserting Eqs. (JTSJ, ipO|). (|2TJ and into Eq. JTTJ), and multiplying by 

(ipm | , yields 

Cn = "^TifTtf ' " = 1.2,3,... (23) 
1 + (k„A_d) 2 



Finally, using Eqs. JTH|), HOJ, © and (23 in Eq. lUHJ) leads to 

a " = fe4 _ ^ 1 + (^) 2 J^< - = 1 ' 2 ' 3 ' 



(24) 



5.^. Transport coefficients 

The flow rate and the electrical current are conveniently written as 

Q = (l\v), (25) 

I ={ P l q \v)+cT E(l\l), (26) 

where the second relation is the linearized Nernst-Planck equation with the first 
term being the convection/streaming current while the second is the ohmic current. 
Substituting Eqs. (|19|) and l|21|) into these expressions we identify the transport 
coefficients as 

A ^ 1 A 

g "-4e^4 l ' ( 27 ) 



n=l n 



G 



<A ^ i At 

ryL ^ 1 + ( Kn A D ) 2 .4 ' 1 ' 

n—l v 7 



12 



G 



21 



where 



- ^ (291 

TyL ^ l + ( Kn A D ) 2 .4 ' 1 j 

n—l v y 

g22 = ^ + ^E KAd)2 A (30) 

A^^^=i<ii^)r of) 

\Wn\Wn) 

is the effective area of the eigenfunction The fraction A n /A is consequently a 

measure of the relative area occupied by \ip n ) satisfying the sum-rule ^™ =1 j4„ = 
A [T]]]. We note that as expected G obeys the Onsager relation G\2 — Gi\- 
Furthermore, using that 

(k„Au) 2 \ d d 1 



[1 + (k„Az,) 2 ] 2 2 d\ D l + {n n \ D ) 2 ' 



(32) 



we get the following bound between the off-diagonal elements G\2 — G21 and the lower 
diagonal element G22, 

a Q A e( dG 12 . . 

G22 --T + 2A7 9A^- (33) 
In the context of the geometrical correction factor a studied in detail in |S] we 
note that the first diagonal element may be written as Gu = G^/a where G\ x = ^ 
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; 2 
k{ 


Ai /A 


a 


In 

a/C 


circle 




4/~/'( ~ 0.69 a ' b 


4tt 


2 C 


quarter-circle 


1.27 d 


0.65 d 


29.97 d 


1.85 d 


half-circle 


1.38 d 


0.64 d 


33.17 d 


l.97 d 


ellipse(l:2) 


1.50 d 


0.67 d 


10tt c 


2.lO d 


ellipse(l:3) 
ellipse(l:4) 


1.54 d 
1.57 d 


0.62 d 
0.58 d 


40tt /3 c 
17tt c 


2.21 d 
2.28 d 


triangle (1:1:1) 


TT 2 /9 ~ 1.10 e 


6/7r 2 ~ 0.61 e 


20\/3 c 


5/3 ~ 1.67 c 


trianp-W 1 -1 -\/2) 


57f2 ~ 1 Dfi a 

4(2+V5)2 - i - U0 


51 2/Q-7T 4 ~ D 58 a 


38.33 d 


1.64 d 


square(l:l) 


7T 2 /8 ~ 1.23 a 


64/tt 4 ~ 0.66° 


28.45 d 


1.78" 


rectangle(l:2) 


5tt 2 /36 ~ 1.37 a 
5tt 2 /32 ~ 1.54 a 
17tt 2 /100~ 1.68 a 
- tt 2 /4 ~ 2.47 a 


64/tt 4 ~ 0.66° 


34.98 d 


1.94 d 


rectangle(l:3) 


64/tt 4 ~ 0.66° 


45.57 d 


2.14 d 


rectangle(l:4) 


64/tt 4 ~ 0.66° 


56.98 d 


2.28 d 


rectangle(l:oo) 


64/tt 4 ~ 0.66° 


oo 


~3' 


pentagon 


1.30 d 


0.67 d 


26.77 d 


1.84 d 


hexagon 


1.34 d 


0.68 d 


26.08" 


1.88" 



Table 1. Central dimcnsionless parameters for different geometries. 
a See e.g. 1241 for the eigenmodes and eigenspectrum. 

'Here, 71 ~ 2.405 is the first root of the zeroth Bessel function of the first kind. 
c See 0. 

d Data obtained by finite-element simulations 1261 . 
e See e.g. |27| for the eigenmodes and eigenspectrum. 
f See e.g. 1221 for a solution of the Poisson equation. 



is a characteristic hydraulic conductance and the geometrical correction factor a can 
be expressed as 

where k n = n n A/V is a dimensionless eigenvalue. In passing we furthermore note 
that this formal result is a convenient starting point for perturbative analysis of the 
correction due to small changes in the boundary dil |25|. 

3.2.1. Non-overlapping, thin Debye layers. For the off-diagonal elements of G we 
use that [1 + {K n X D ) 2 ]~ l = 1 + £>[fc 2 (X D V/A) 2 ]. In Section El we numerically 
justify that the smallest dimensionless eigenvalue k\ is of the order unity, so we may 
approximate the sum by a factor of unity, see Tabled whereby we arrive at Eq. JSJ 
for Xd <C AjV . These results for the off-diagonal elements are fully equivalent to the 
Helmholtz-SmoluchowsM result [IS] . For G22 we use that (k„Az}) 2 [1 + (K n Ai)) 2 ]~ 2 = 
C[/c 2 (XdV/A) 2 ], thus we may neglect the second term, whereby we arrive at Eq. J2J|. 

3.2.2. Strongly overlapping Debye layers. In the case of k,\Xd > 1 we may use the 
result [l + {KnXo) 2 }- 1 = {n n X D )- 2 + 0[k-\X D V /A)- 4 } which gives 

eC A 
G12 = G21 — — -rj- Gn, Xd ^> — ■ (35) 
A D F 
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This is the Debye-Hiickel limit of Eq. (J5J for strongly overlapping Debye layers. For 
G22 we use Eq. (|33H and arrive at the result in Eq. © for ZeQ <C k B T by using 



eD/X 2 



D ' 



3.2.3. The circular case. For a circular cross-section, where a = 87r, it can be shown 
that fig) 



(36) 



i V h^SA/aX 2 

where /„ is the nth modified Bessel function of the first kind, and were we have 
explicitly introduced the variable A/a\ 2 D to emphasize the asymptotic dependence in 
Eq. for strongly overlapping Debye layers. We note that we recover the limits in 
Eqs. © and @ for A fl « A/V and \ D > A/V, respectively. 



4. Beyond the Debye Hiickel approximation 

In order to go beyond the Debye-Hiickel approximation we consider, for simplicity, a 
symmetric binary (Z:Z) electrolyte. Next, we neglect strong correlations between the 
ions so that the equilibrium properties of the electrolyte are governed by Boltzmann 
statistics [H], i.e., the concentrations of the two type of ions are given by 

Ze 



c eq(n-) = c Q exp 



k R T 



<l>eq(r±) 



(37) 



This is equivalent to assuming equilibrium with bulk reservoirs at the ends of the 
channel in which the potential (j) eq tends to zero and both concentrations c~ q to c D . 
Substituting the charge density pi = Ze(c+ —c~) into the Poisson equation ((TTIl 



we arrive at the nonlinear Poisson-Boltzmann equation [31 118| 



_LVeq 



(rx) - 



k B T 
Ze\ 2 D 



sinh 



Ze 



^eq(rj_) 



where the Debye screening length is given by 



ek B T 



V 2{Ze) 2 c a 

The conductivity er of the electrolyte depends on the local ionic concentrations 



(Ze) 2 D 
k B T 



"cq 



(rx) 



"cq 



(fj_)] = (7 COsh 



Ze 
k^f 



^eq(rL) 



(38) 



(39) 



(40) 



assuming equal diffusivities D for the two ionic species. In the Debye-Hiickel limit, 
Ze£ <C k B T, where thermal energy dominates over the electrostatic energy we may 
linearize the right-hand side of Eq. (|38|) so that we arrive at the Debye-Hiickel 
approximation in Eq. |JTBJ|. Similarly the conductivity in Eq. I|4U|I reduces to the 
bulk conductivity er . Compared to Eq. (|26|l the electrical current obtained from the 
nonlinear Nernst-Planck equation becomes 

I=(pM+E(a\l). (41) 

We calculate the off-diagonal elements from G12 = G21 = (pl q \v p ) / Ap and find 

G12 = G21 



ek B T I 1 Ze 
— ^-( smn 

ApZe\ 2 D \ \k B T 



(42) 
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Similarly, Eq. I|15|) for the two components in the electrical conduction G22 

G"^ + G conv we get 

(Jo I , / Ze 

= l\ cos H 



sin 

(JTO 



1 , 



Gconv 
22 



2L k B T 



k B T 



(43) 
(44) 



where we have used that <7 = eD/\ 2 D and introduced the dimensionless quantity m, 



(k B T\ 2 2e 
\ Ze J r}D' 

which indicates the importance of electro-osmosis relative to electro-migration. 



(45) 



4-1. Non- overlapping, thin Debye layers 

In the limit of thin Debye layers we have already discussed how Eq. (|14fl in general 
leads to Eq. @ because the screening is good and (f> cq is nonzero only on a negligible 
part of f2. This property is more implicit when G12 or G21 is written in the form 
of Eq. (|42|l . which is a more appropriate starting point for analyzing the limit of 
strongly overlapping Debye layers. For G22 the calculations are more involved; we 
assume that the channel wall is sufficiently smooth on the Debye-length scale so that 
we can everywhere use the Gouy-Chapman (GC) solution for a semi-infinite planar 
geometry [31 H%]. 

x gc _ k BT_ 

where r n denotes the normal distance to the channel wall. Substituting this into 
Eq. (|43|) and l|44|) the integrals can be carried out analytically resulting in 

G22 = *>j(l + Du), A D «^ ; (47) 
where Du is the Dukhin number 

^ = ^(l + m)4sinh 2 f-f^Y (48) 



Unh (l?4=) e -m/Xn 



(46) 



A v ; \4:k B T 

defined as the ratio of the surface conductivity in the charged Debye layers to the bulk 
conductivity a times the geometrical length scale A/V (see Ref. |2E1 and references 
therein). Clearly, when the Debye layer becomes very thin, surface conduction is 
negligible and we recover the simple result in Eq. Q. 



4-2. Strongly overlapping Debye layers 

When the Debye layers are strongly overlapping the screening is weak and 4> cq w ( 
throughout the cross section. Hence we can pull the integrand smh(Ze(j) cq /k B T) 
outside the bra-ket in Eq. I|42fl and we arrive at Eq. (J5J. Here, we have used that 
(l|u p )/Ap = G11 and introduced the parameterization in Eq. J2J. Similarly, from 
Eqs. I|43|l and H44J) we obtain Eq. © where we have used Eqs. (|14|l and Q to eliminate 

(i|C 

We note that due to shifts in free energies, the zeta potential inside a narrow 
channel with significant Debye-layer overlap is generally not the same as in a 
macroscopic channel with no overlap, see e.g. ^H2U for a discussion. 
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5. Numerical results 

5.1. The Helmholtz basis 

Only few geometries allow analytical solutions of both the Helmholtz equation and 
the Poisson equation. The circle is of course among the most well-known solutions 
and the equilateral triangle is another example. However, in general the equations 
have to be solved numerically, and for this purpose we have used the commercially 
available finite-element software Femlab j^HJ. The first eigenstate of the Helmholtz 
equation is in general non-degenerate and numbers for a selection of geometries are 
tabulated in Table 1. Note how the different numbers converge when going through the 
regular polygons starting from the equilateral triangle through the square, the regular 
pentagon, and the regular hexagon to the circle. In general, k\ is of the order unity, 
and for relevant high-order modes (those with a nonzero An) the eigenvalue is typically 
much larger. Similarly, for the effective area we find that Ai/A < 4/ 7l 2 ~ 0.69 and 
consequently we have A n /A < 1 — 4/7^ ~ 0.31 for n > 2. The transport coefficients 
in Eqs. i|27|l to l|30|l are thus strongly influenced by the first eigenmode which may be 
used for approximations and estimates of the transport coefficients. As an example 
the column for a/C is well approximated by only including the first eigenvalue in the 
summation in Eq. I|34|) . 



5.2. Transport coefficients 

Our analytical results predict that when going to either of the limits of thin non- 
overlapping or strongly overlapping Debye layers, the transport coefficients only 




Figure 2. Rcscaled off-diagonal transport coefficients versus rescaled Debye-layer 
thickness in the Debye-Hiickel limit. The solid line is the exact result for a circle, 
Eq. 1361 . and the dashed line shows Eq. JSJ for Zef -C fcsT. The data points are 
finite-element simulations in the linearized regime for different cross sections, see 
inset. Finally, the dash-dotted line shows the first term from the summation in 
Eq. only. 
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Figure 3. Rescaled off-diagonal transport coefficients versus rescaled Debye- 
layer thickness beyond the Debye— Hiickcl approximation. The solid line is the 
exact result for a circle within the Debye— Hiickel approximation, Eq. 1361 . The 
data points are finite-element simulations for different cross sections (see inset) 
with ZeC^/kgT = 1, 3, 6, and 10 from below. The dashed lines indicate the 
corresponding asymptotic expression for strong Debye-layer overlap, Eq. l5l . 



depend on the channel geometry through the cross sectional area A and the correction 
factor a. Therefore, when plotted against the rescaled Debye length \Jctj A Xd, all 
our results should collapse on the same asymptotes in the two limits. 

In Fig. we show the results for the off-diagonal coefficients obtained from finite- 
element simulations in the Debye-Hiickel limit ZeQ <C fcsT for three different channel 
cross sections, namely two parabola shaped channels of aspect ratio 1:1 and 1:5, 
respectively, and a rectangular channel of aspect ratio 1:5. In all cases we find excellent 
agreement between the numerics and the asymptotic expressions. For the comparison 
we have also included exact results, Eq. (J2EJ), for the circular cross section as well as 
results based on only the first eigenvalue in Eq. (|2EJ). Even though Eq. is derived 
for a circular geometry we find that it also accounts remarkably well for even highly 
non-circular geometries in the intermediate regime of weakly overlapping Debye layers. 

In Fig. |3 we show numerical results for the off-diagonal transport coefficients 
beyond the Debye-Hiickel approximation. At large zeta potentials the Debye layer 
is strongly compressed and the effective screening length reduced. Therefore the 
suppression of the electro-osmotic flow/streaming current at strong Debye-layer 
overlap is shifted to larger values of Ad as compared to the Debye-Hiickel limit. For 
the comparison we have also included the exact result for a circular cross section in the 
Debye-Hiickel approximation as well as the asymptotic expression for non-overlapping 
and strongly overlapping Debye layers, Eq. As seen the asymptotic expressions 
account well for the full numerical solutions independently of the geometry. 

In Fig. 01 we show numerical results for the electrical conductance beyond the 
Debye-Hiickel approximation. Open symbols show the electro- migration part G^' s , 
subtracted the trivial bulk contribution a A/L, whereas solid symbols show the 
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Figure 4. Comparison of the two components in the G22 transport coefficient 
for three different geometries, see inset. Open symbols show (G^ ls L/cr ,4) — f, 
i.e., the surface specific contribution to the electrical conductance from electro- 
migration, see Eq. 1431 . Solid symbols show the contribution from electro-osmotic 
convection G^ 1 ^ L / ma A, see Eq. 1441 : for ease of comparison we have included 
the dimensionless number m in the rescaling. The data are obtained from finite- 
element simulations for different cross sections (see inset) with Ze^/ksT = 1, 3, 
6, and fO. The solid lines indicate the corresponding asymptotic expressions 
for strong and weak Debye-layer overlap, Eqs. JgJ and 1471 . respectively, for 
Ze^/kgT = f in a circular channel. Asymptotes for the other values of ZeC^/kgT 
have been left out for clarify, but equally good agreement is found also in these 
cases. 



electro-osmotic convection part , see Eqs. (|43|) and (|44|) . respectively. Again, 

we find that the numerics are in excellent agreement with our asymptotic results. For 
the Xd <^V/A regime we note that the Dukhin number, Eq. is proportional to 
XdV/A = \/Cj 'A\d and not yjaj AXu- Therefore, strictly one would not in general 
expect data to collapse on the same asymptote for Xd <C V j A. Looking carefully at 
this part of the figure, small variations can be seen from geometry to geometry - the 
reason why the variations are still so small is that a/C ~ 2 independently of geometry, 
see Table 1. 



6. Conclusion 



We have analyzed the flow of incompressible electrolytes in long, straight channels 
driven by pressure and electro-osmosis. By using a powerful Hilbert space 
eigenfunction expansion we have been able to address the general problem of an 
arbitrary cross section and obtained general results for the hydraulic and electrical 
transport coefficients. Results for strongly overlapping and thin, non-overlapping 
Debye layers are particular simple, and from these analytical results we have calculated 
the corrections to the hydraulic resistance due to electro- hydrodynamic interactions. 
These analytical results reveal that the geometry dependence only appears through 
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the area A and the correction factor a, as the expressions only depend on the rescaled 
Debye length ^/a/A Xd- Our numerical analysis based on finite-element simulations 
indicates that these conclusions are generally valid also for intermediate values of 
Arj. Combined with recent detailed work on the geometrical correction factor [S] the 
present results constitute an important step toward circuit analysis ^5] of complicated 
micro and nanofluidic networks incorporating complicated cross-sectional channel 
geometries. 
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